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We study the sensitivity of the higher-order moments of produced particle multiplicity distri¬ 
butions to the chemical freeze-out parameters in relativistic heavy ion collisions using the Hadron 
Resonance Gas (HRG) model. We compare the obtained sensitivity level to the one extracted from 
the ratios of particle yields. We find that, for certain final state hadrons, the fluctuation measure¬ 
ments add significant information to the determination of the hadro-chemical freeze-out properties 
of the deconfined phase of matter obtained at RHIG and the LHG. 


INTRODUCTION 

Over the past few decades the particle production in 
relativistic heavy-ion collisions has been successfully de¬ 
scribed by statistical hadronization models which postu¬ 
late a chemical freeze-out surface common to all hadronic 
species. This means the yield of particles emitted from 
the interaction region can be reproduced by assuming an 
equilibrated system with fixed thermodynamic param¬ 
eters. At this time in the dynamical evolution of the 
expanding fireball, the production of particles through 
inelastic collisions ceases and freeze-out of the hadro- 
chemistry is achieved. The resulting phase space surface 
is generally characterized (at sufficiently high collision 
energies) by only two free parameters, namely the che¬ 
mical freeze-out temperature and the associated baryo- 
chemical potential {Tch and TB,ch)- These types of static 
models have proven to be extremely successful over a 
wide range of collision energies, which led to the defini¬ 
tion of a generally accepted chemical freeze-out curve in 
the T — piB plane [1] . It was also shown that this curve 
is close to the QCD transition curve as defined by lat¬ 
tice QCD [5]. Since the QCD transition is an analytic 
crossover for small pis [3], the lattice curve corresponds 
to the steepest gradient in the temperature dependence 
of some characteristic observables which identify the de¬ 
confinement and chiral restoration transition nil]. 

In the past, statistical model calculations were used 
only to describe the multiplicity of particles, either by 
calculating the expected yields (which adds the interac¬ 
tion volume as an additional free parameter to the cal¬ 
culation) or by calculating particle ratios (in which case 
the volume cancels out, assuming a common, flavor in¬ 
dependent freeze-out volume) [TJ |3H3| . 

Recently, it has been shown that the moments of net- 
particle multiplicity distributions from the experiment 
can be related to susceptibilities of conserved charges 
calculated on the lattice unHia. This allows the di¬ 
rect determination of chemical freeze-out parameters in 
the thermally equilibrated grand-canonical ensemble ap¬ 
proach on the lattice without having to rely on statistical 
models. 


The continued application of statistical models in the 
context of these lattice QCD-to-experiment comparisons, 
is based on the fact that conserved charges on the lattice 
can only be directly related to particle distribution mo¬ 
ment measurements by studying and controlling the lim¬ 
itations of the experiment to measure conserved charges. 
Thus, by modeling detector effects, such as acceptance 
and reconstruction efficiency, in a statistical model cal¬ 
culation, the latter de-facto facilitates the necessary link 
between experiment and lattice QCD. It was shown that 
statistical Hadron Resonance Gas (HRG) models repro¬ 
duce the equilibrium lattice QCD results for the low¬ 
est order susceptibilities and their ratios in the hadronic 
phase reasonably well |I2HT4] . An added advantage of 
the HRG model calculations is the fact that a conserved 
charge, such as strangeness, can be broken down into its 
contributions from each identified particle species, which 
is not achievable on the lattice. 

We have shown previously that the calculations of the 
lowest moments of particle multiplicity distributions in 
a statistical model are in agreement with recent STAR 
data [15]. One can therefore ask whether the chemi¬ 
cal freeze-out parameters obtained from our analysis are 
comparable to parameters from an analysis of the particle 
multiplicities. For this purpose, we study here the sen¬ 
sitivity of higher-order moment ratios on the freeze-out 
parameters in comparison to that of multiplicity ratios. 
Depending on the observable and the experimentally 
achievable accuracy in the corresponding measurement, 
it might be advantageous to consider fluctuations rather 
than yields for pinning down the chemical freeze-out con¬ 
ditions. This requires that the measured fluctuations are 
those of an equilibrated hadronic medium [Min]- 
From the experimental data the central moments of 
the (net-)particle distributions can be constructed. The 
corresponding cumulants are given by: 

Cl = ((7V)> C2 = {{5Nf) 

C3 = {{6Nf) C4 = {{SNf) - 3{{5Nff 

where SN = N — (N) is the fluctuation of the (net- 
)particle multiplicity around its mean. The cumulants Ci 
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relate directly to the susceptibilities Xii which are quan¬ 
tities that can be calculated for thermodynamic systems, 
e.g. in lattice QCD. Susceptibilities are defined as deriva¬ 
tives of the pressure with respect to the chemical poten¬ 
tial. In order to quantify certain features of distributions 
beyond the mean (M = xi) and the variance (cr^ = X 2 ) 
one often looks at the skewness S and the kurtosis k de¬ 
fined as: 

mean : M = xi variance : a-^ = X 2 
skewness : S = X'ilX'J'^ kurtosis : k = X 4/X2 

Then, one can relate the thermodynamic susceptibi¬ 
lities calculated on the lattice or in an HRG model to 
quantities obtained experimentally through the following 
volume-independent ratios: 

S<T = X3lX2 = X 4 /X 2 

Mla'^ = XllX2 Sa^/M = X3/Xl- 


DETAILS OF THE HADRON RESONANCE GAS 
MODEL 

Detailed measurements of conserved charge fluctua¬ 
tions have been conducted at several collision energies 
at RHIC in order to search for non-statistical fluctua¬ 
tions that could signal the existence of a critical point 
in the QCD phase diagram [T8ll20j . A series of studies 
using variations of the standard statistical hadronization 
models to determine the baseline for these critical point 
searches can be found in the literature The re¬ 

sults presented in this paper are obtained using an HRG 
model in partial chemical equilibrium, which means that 
all contributions from the strong decay of hadronic res¬ 
onances are taken into account. The list of all included 
resonant states is based on the PDG tables [26], up to a 
mass of 2 GeV/c^. 

Further details of this HRG model can be found in m, 
where our group compared the net-charge and net-proton 
distribution measurements from the STAR collaboration 
to the HRG model results. We evaluated ratios of the 
lower moments of these multiplicity distributions as func¬ 
tions of the temperature T, baryo-chemical potential 
and the chemical potentials /rg and /X 5 . The relation 
among these thermodynamic variables is obtained by im¬ 
posing certain initial conditions occurring at the colli¬ 
sion, namely the conservation of the total net-strangeness 
ns = 0 and the proper ratio of protons to baryons in the 
colliding nuclei ng = ^ub- 

We obtained lower freeze-out temperatures (wldG 
MeV) than the statistical hadronization model fits to 


all available particle multiplicities («166 MeV). We at¬ 
tributed this discrepancy in part to the fact that net- 
proton and net-charge are mainly sensitive to the light 
quark susceptibilities, whereas a full fit to all parti¬ 
cle multiplicities contains a significant contribution from 
(multi-)strange particles. A possible separation of light 
and strange quark transitions was suggested by high pre¬ 
cision lattice QGD simulations of the kurtosis [27]. If 
true, then strange particle multiplicities and moments 
might yield a higher freeze-out temperature, as already 
pointed out in [TS] HH] ■ 

A comparison of our freeze-out temperature obtained 
from the net-proton and net-charge distribution to the 
one from measured particle yields (see Ref. [H]), indeed 
indicates a discrepancy between the temperature neces¬ 
sary to describe the strange and multi-strange baryon 
yields compared to the lower mass and mostly light quark 
states (tt, AT, p). In this fit to the particle multiplicities 
pions and kaons, though, show very little sensitivity to 
the chosen freeze-out temperature at fixed chemical po¬ 
tential. Since pions, kaons and protons are presently the 
only particles for which the experiments can determine 
higher order fluctuations, it is interesting to ask whether, 
for these hadronic final states, the moments of the net dis¬ 
tributions are more sensitive to the chemical freeze-out 
conditions than the basic particle ratios. 

Our initial study has also shown that the highest cu- 
mulants might be prone to non-statistical effects such as 
volume fluctuations or chiral criticality, even at collision 
energies in the crossover region |3D|. Critical behavior 
might be captured in lattice QCD, but statistical models 
are by definition insensitive to these dynamical fluctua¬ 
tions. Therefore we initially focus our study on the lowest 
moments necessary to determine freeze-out parameters, 
namely X 2 and xi ■ Whether ^3 already contains a small 
contribution from critical fluctuations is an open issue, 
which is presently debated in the literature. For very high 
collision energies, well separated from a potential critical 
point, these effects should be small. Thus, for the highest 
RHIC energies and all LHC energies it should be safe to 
also use xa or even X 4 for a potential determination of 
the chemical freeze-out parameters, and we show results 
also for moment ratios that include the higher moments. 


RESULTS 

Based on the above mentioned tension in the extracted 
chemical freeze-out temperatures using the STAR net- 
proton and net-charge fluctuation results compared to 
the particle ratios, in particular at small baryo-chemical 
potentials |15j . we perform the comparative study pre¬ 
sented here only for the highest RHIC energy. Since 
there is good agreement between the baryo-chemical po¬ 
tentials deduced from the particle ratios and the fluc¬ 
tuations, we fixed at 24.3 MeV and determined the 
sensitivity of different observables to the freeze-out tem¬ 
perature. Historically, ratios of mesons and baryons to 
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pions have proven to be a good thermometer, since their 
dependence on the temperature T^h is mainly driven by 
the particle mass difference. Particle yields might have a 
higher sensitivity to the temperature but a fit to ratios is 
preferred, since the volume cancels out and the ratios are 
less prone to biases [3T] . Our results for the sensitivity of 
particle ratios to the temperature are in good agreement 
with previous statistical model calculations [3TJ[35]. Here 
we add the sensitivity of the moments to this study, and 
we show that for certain particles it might be beneficial 
to use the moment analysis, if efficiency corrected data 
are available. 

A comparison of the freeze-out temperature sensitiv¬ 
ity of the particle ratio (with respect to the pion yield) 
and the lowest moment ratio (X 2 /X 1 ) for all main hadron 
species {K,p,A, S, fl) is shown in Fig. In particular 
for kaons the difference is striking: the AT/tt ratio is al¬ 
most completely insensitive to the temperature, whereas 
the net-kaon moment ratio shows a strong sensitivity well 
outside the achievable experimental error bars. 

The experimental error bars, shown as vertical thick 
lines in Fig. 1, are from for the particle ratios, the 
net-proton X 2 /X 1 uncertainty is from the recent STAR 
publication [TB]. For the net-Kaon X 2 /X 1 error we have 
used a preliminary study by STAR [5D] based on mo¬ 
ments that were not efficiency corrected. The effect of 
the reconstruction efficiency correction was estimated to 
roughly double the uncertainty |35j and the error bar 
shown in Fig. 1 reflects the anticipated error for an ef¬ 
ficiency corrected measurement. In general, we show 
experimental uncertainties only for particle species for 
which both measurements (the yield ratios and the net- 
moment ratios) have been published. 

It is interesting to note that for protons, which re¬ 
quire a lower temperature for the yields and moments at 
the LHC and the higher RHIC energies, both measure¬ 
ments show the necessary sensitivity based on the exper¬ 
imentally achievable error bars. The lower temperature 
deduced for protons from lattice QCD and HRG model 
moment analyses nmn], might thus explain the ’proton 
anomaly’ in the yield. Generally, the baryon yield ra¬ 
tios show a stronger sensitivity to the temperature than 
the corresponding lower moment ratios; nevertheless, for 
multi-strange hyperons, both particle yield and moment 
ratio exhibit a comparable steepness in the relevant tem¬ 
perature range. 

The moment ratios which we show in Fig. 1 include 
kinematic cuts in rapidity and transverse momentum ac¬ 
cording to the experimental restrictions as well as contri¬ 
butions from resonance decays. For the latter, two effects 
have to be taken into account. First, a resonance can 
decay and be regenerated, which leads, for sufficiently 
short-lived resonances, to an isospin randomization for 
the daughter hadrons. This effect is for example relevant 
for the net-protons since the A-resonances have a very 
short lifetime. Kitazawa and Asakawa [36l |37j have pro¬ 
posed corrections, which are applied to the net-proton 
moment ratio in Fig. 1 (see also [15]). The magnitude of 
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FIG. 1. (Color online) Comparison of the sensitivity of Xz/Xi 
ratios and particle yield ratios to the chemical freeze-out tem¬ 
perature. All values are normalized to the value at 140 MeV. 
The baryo-chemical potential has been fixed to 24.3 MeV. The 
vertical bars in the plots show the presently achievable exper¬ 
imental uncertainties for particle species for which both par¬ 
ticle yield and moment ratios are available (based on STAR 
measurements [l8l [20l [33l [3^ 1. 


the kinematic cuts and the isospin randomization effect 
on the slope of X 2 /X 1 is shown in Fig. 2 for net-protons: 
while different cuts on momentum and rapidity have a 
very small effect [23], the corrections due to isospin ran¬ 
domization (KA corrections) visibly increase the sensi¬ 
tivity of this curve to the temperature. The net-proton 
X 2 /X 1 shown in Fig. 1 corresponds to the red, dotted 
line in Fig. 2, for which all possible effects are taken into 
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FIG. 2. (Color online) Net-proton X 2 /X 1 a function of the 
temperature: comparison between the curves obtained with 
different kinematic cuts and with/without isospin random¬ 
ization. All values are normalized to the proper value at 140 
MeV. 


account. 

The second relevant effect is due to the probabilistic 
nature of the resonance decay, because the number of 
ground state hadrons follows a probability distribution 
due to the different branching options. This was first 
pointed out in m and can be applied to the suscepti¬ 
bilities Xi with I > 1 by taking into account the actual 
decay contributions rather than an averaged branching 
pattern [23]. We apply this probabilistic effect to all par¬ 
ticles investigated in Fig. 1 except the (anti-)protons as 
here isospin randomization dominates the probabilistic 
resonance decay contributions. Fig. shows the net- 
kaon and net-A X 2 /X 1 with and without the probabilis¬ 
tic decay contribution (PROB). As in the case of the 
isospin corrections, the temperature sensitivity changes 
when the effect is taken into account. It is interesting 
to note, though, that the relative temperature sensitivity 
for net-kaon fluctuations reduces, whereas for net-As the 
sensitivity is enhanced. 

The upper panel of Figj^ shows a comparison of the 
net-X 2 /xi temperature sensitivity for all relevant stable 
hadrons that are measurable by the experiments. A no¬ 
ticeable baryon-to-meson difference can be established. 
This difference can be understood qualitatively already 
for primordial hadrons (without resonance decays) in the 
Boltzmann approximation, where the net-cumulant ratio 
X 2 /X 1 is given by coth(/ii/T), with being the che¬ 
mical potential of particle i jUj. The ratio ^ has a 
different trend for mesons and baryons, as can be seen 
in the lower panel of Fig. (note that the normalized 
difference from the value at T = 140 MeV is shown), be¬ 
cause the individual components that contribute to /ii, 
i.e. pLB, piQ, sum up to an opposite temperature de¬ 
pendence for baryons and mesons. 

Let us point out that, due to resonances which can 
decay into a meson and/or its anti-meson, the suscep¬ 
tibilities for net-mesons cannot simply be obtained via 
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FIG. 3. (Golor online) Net-kaon and net-A Xa/xi s-s S' 
function of the temperature: comparison between the curves 
with/without the probabilistic decay contribution (PROB). 
All values are normalized to the proper value at 140 MeV. 


independent production. For the second susceptibility 
this is seen from: 

X2,net-i =X2 + X2* + 

- 2{n,)R{ni)R)x2^- ( 1 ) 

R 

The first two terms, X 2 *i represent the contributions 
from the primordial particles, while the sum in Eq. ([^ 
accounts for the average resonance decay contributions. 
The last term in Eq. ([^ is non-zero whenever a res¬ 
onance decays into a stable particle, its antiparticle or 
the particle-antiparticle pair. Since there are no known 
hadronic resonances that decay into a baryon and/or its 
anti-baryon, Eq. ([^ reduces to independent production 
in the case of baryons, while correlations arise for mesons. 

It has been established, on the basis of simple parti¬ 
cle yields, that in a statistical model the strange baryons 
are more sensitive to the temperature since the strange 
quark rest mass is close to the temperature of the equi¬ 
librated system. As shown in Eig. [^ the steepness of 
the X 2 /X 1 curves, and thus the sensitivity to the temper¬ 
ature, increases with the strange quark content, which 
means that cumulant ratios in the multi-strange baryon 
sector rise significantly faster than for the A. This effect 
is mainly caused by imposing vanishing net-strangeness, 
which means that even at constant the strange che¬ 
mical potential fj,s has to increase as a function of the 
temperature (in our simulation from 3 to 8 MeV between 
140 and 180 MeV). The impact can already be seen in 
the temperature dependence of the A X 2 /xi- For the 
multi-strange baryons this rise grows according to the 
strange quark content, each quark leading to an addi¬ 
tive strangeness chemical potential contribution. In other 
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FIG. 4. (Color online) Upper panel: comparison of the sensi¬ 
tivity of ^ ratios to the chemical freeze-out temperature for 

different hadron species. Lower panel: for particle 

i as a function of the temperature for mesons and baryons. 
All values are normalized to the proper value at 140 MeV and 
the baryo-chemical potential has been fixed to 24.3 MeV. 


words, for the net-fl at high T the strangeness chemi¬ 
cal potential (/rs = 3x8 MeV = 24 MeV), offsets the 
baryo-chemical potential {fJ-B = 24.3 MeV) since the two 
chemical potentials always carry opposite sign. At this 
point the net-multiplicity grows much slower than the 
variance, which leads to a strongly increasing X 2 /X 1 
tio. This effect is not unique to the specific value of 
the baryo-chemical potential, in other words at lower 
and higher collision energies the temperature dependence 
of the strange baryons will have similar trends to the 
ones obtained here. The relative increase of Hs{T) in 
our HRG model is comparable to the one deduced from 
lattice QCD, but the absolute values do not match the 
dependence from lattice QCD [T^], which is about 15% 
higher than in a standard HRG model. We note that 
the splitting between the net-pion and net-kaon results 
is, likewise, a consequence of the different strangeness 
content and the imposed strangeness neutrality. 

Preliminary studies, using all measured resonance 
states listed in PDG-2014 [3S] up to a mass of 2.5 
GeV/c^, show that the inclusion of these known higher 
mass resonances further improves the agreement with 
the lattice QCD [39]. This might suggest the 

existence of even higher mass, as of yet unidentified, 
strange states, some of which are predicted by quark 
models gni mi, but any further expansion of the res¬ 
onance mass spectrum needs to be rigorously tested 
against all sensitive susceptibility predictions from lattice 
QCD. We have not included those Hagedorn type states 
in our calculation. We have verified, though, that the 
higher fis caused by the additional states listed in PDG- 
2014 will further increase the temperature dependence 
of the multi-strange baryon X 2 /X 1 shown in Fig. El(by 


(10 — 25)% in the relevant temperature range), whereas 
the particle yield ratios stay rather unaffected. 
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FIG. 5. (Golor online) Comparison of sensitivity of the net- 
Xa/Xi, X 3 /X 2 and X 4 /X 2 ratios to the chemical freeze-out 
temperature. All values are normalized to the proper value 
at 140 MeV. The baryo-chemical potential has been fixed to 
24.3 MeV. 

Fig. 0 shows that the sensitivity of the lower cumu- 
lant ratios is also preserved in the X'ilxi ratio, although 
the sign of the temperature dependence changes. As ex¬ 
pected, the kurtosis-based ratio (X 4 /X 2 ) shows signifi¬ 
cantly less sensitivity, but its determination is crucial to 
obtain the strange quark freeze-out conditions, since in 
lattice QCD, only the even strangeness susceptibilities 
yield a finite result. Therefore any link to first princi¬ 
ple calculations needs to be established through X 4 /X 2 , 
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which is the lowest even ratio that can be measured. As 
shown in Fig. each strange particle by itself shows 
very little temperature dependence, but the more sensi¬ 
tive total net-strangeness X4 X 2 can be calculated 
using a combination of all contributing strange particles 
with pre-factors based on their strangeness content. It 
should be noted that no PROB corrections were applied 
to the net-higher moment ratios in Figj^ 



FIG. 6. (Color online) Upper panel: comparison of 
^et-s j^net-s obtained from different strange particle com¬ 
binations. The black squares are the lattice QCD results from 
Ref. |27) . Lower panel: comparison of ob¬ 

tained from different strange particle combinations. 

In Figj^we address the question whether measuring a 
subset of particle species that carry strangeness is suffi¬ 
cient to relate the experimental results to calculations of 
the full strange quark content on the lattice. In the upper 
panel, we show a comparison of Xa ^*'~^obtained 
from different strange particle combinations in our HRG 
model. It is evident that, in order to achieve a sensitivity 
to the temperature, which is comparable to the one on 
the lattice, multi-strange baryons need to be included in 
the analysis. Obtaining sufficient statistics in the exper¬ 
iment for higher moments becomes increasingly difficult 
for such baryons. 


A study purely based on a comparison between exper¬ 
imental data and HRG model predictions can be based 
on X 2 /A 1 alone, though, since there is no advantage in 
using the higher cumulants in terms of sensitivity. The 
lower panel of Fig. shows that, contrary to our higher 
moment study for X 4 /X 2 , the sensitivity of X 2 /X 1 does 
not benefit from the inclusion of multi-strange baryons: 
all curves are parallel, thus yielding a comparable sen¬ 
sitivity to T. Given the difficulty of obtaining sufficient 
statistics for multi-strange baryons, it might therefore be 
even advantageous to look only at lighter strange hadrons 
as long as one only wants to compare to statistical model 
predictions. 

CONCLUSIONS 

In this paper we have shown that, for certain final state 
hadrons, the statistical analysis of the higher moments of 
the net-particle distributions might add significant sen¬ 
sitivity to the determination of chemical freeze-out pa¬ 
rameters compared to a model that uses particle yields 
and ratios only. In particular, we found that the tem¬ 
perature sensitivity of the net-kaon X 2 /X 1 ratio in com¬ 
parison with the experimentally achievable accuracy is 
sufficient to reliably determine Teh, quite in contrast to 
the K/tt ratio. For protons we found a comparable tem¬ 
perature sensitivity to experimental accuracy relation for 
both, particle and moment ratios. 

Our study indicates that it should be possible to obtain 
a temperature for the chemical freeze-out of strangeness 
from the efficiency corrected net-kaon X2/xi- This might 
be of particular interest if the strange hadron sector dis¬ 
plays a higher Teh than the one extracted from net-proton 
and net-charge distributions, which are dominated by 
light quark particles. For a comprehensive study of the 
chemical freeze-out properties, we suggest that the exper¬ 
iments provide detailed moment analyses for all identified 
particle species with sufficient event-by-event yields. Ef¬ 
ficiency corrected experimental results for the lowest four 
moments are already available for net-protons from the 
RHIG beam energy scan. Uncorrected data have also 
been shown for net-kaons, and we are anticipating cor¬ 
rected net-kaon moment distributions from RHIG and 
LHG in the near future. 
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